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Abstract. Investigating the canonical problem of a periodic multilayered stack 
containing isotropic chiral layers, we homogenized it as a uniaxial bianisotropic medium 
and derived its effective constitutive parameters. The stack shows a resonant behavior, 
when its unit cell consists of a metallic layer and an isotropic chiral layer. The 
presence of isotropic chirality can result in small shifts of the resonance frequency 
for reasonably large values of the chirality parameter, implying that the sign of an 
effective permittivity can be switched. Such spectral shifts in the dielectric properties 
can be potentially useful for spectroscopic purposes. 
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1. Introduction 

Metamaterials that refract negatively [1, 2] have been demonstrated across the 
electromagnetic spectrum from microwave to optical frequencies. A metamaterial's 
response characteristics and frequency range of operation are principally determined 
by the structures of its constituent units. There is, however, a need to be able to 
dynamically tune the response of a given metamaterial once it has been fabricated. 
Kerr nonlinearity [3] and the electro-optic effect [4] have already been suggested in this 
regard. These effects, however, are primarily dielectric effects. A related question arises: 
can isotropic chirality can be utilized to manipulate the properties of a metamaterial 
wherein an isotropic chiral material has been embedded? 

The rotation of the vibration ellipse of light after passage through an isotropic 
chiral material has been known for about two centuries [5]. In the frequency domain, 
this kind of material is described by three constitutive parameters: relative permittivity, 
relative permeability, and chirality parameter. At least one of the three constitutive 
parameters has to be different from its free-space counterpart in order that the material 
be distinguishable from free space (i.e., vacuum). For an isotropic chiral material to 
be negatively refracting [6]- [8], a minimum of two constitutive parameters must be 
different from their free-space counterparts [9, 10]. Thus, the inclusion of isotropic 
chirality enlarges the potential for fabricating isotropic, homogeneous, and negatively 
refracting materials that are actually useful. 

For an isotropic chiral material to be negatively refracting, the chirality parameter 
must be sufficiently large in magnitude compared to the product of the relative 
permittivity and the relative permeability, the underlying assumption being that the 
material is nondissipative at the frequency of interest [8]. If dissipation is present, 
the condition for negative refraction becomes more complicated, but, in essence, the 
magnitude of the chirality parameter is still required to be large [10]. In nature, large 
chirality parameters (in the optical regime) are not known to exist [11], and very large 
chirality parameters for artificial materials (in the microwave regime) have not been 
reported yet [12, 13, 14]. Therefore, for isotropic chirality to be effective in delivering the 
attribute of negative refraction, the magnitude of either the relative permittivity or the 
relative permeability must be close to zero [8]. For dynamic control of a metamaterial 
using isotropic chirality, the fact that the chirality parameter can dominate relative 
permittivity and/or relative permeabilty becomes particularly relevant in the frequency 
ranges where either of them goes through a zero. 

In this communication, we investigate whether the isotropic chirality of an 
embedded medium can be effectively used to strongly influence the spectral properties of 
a metamaterial by choosing the canonical problem of a periodic multilayered stack that 
can be homogenized. Provided the periodic multilayered stack consists of metallic layers 
and isotropic chiral layers, we show that the latter type of layers can even switch the sign 
of the real part of the effective permittivity of the stack. The associated frequency shifts 
are of the order of an angstrom in the optical regime, which are easily measurable [15] 
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and could be useful for spectroscopic purposes. 

The plan of this communication is as follows. Section 2 describes two different 
formulations to homogenize a periodic multilayered stack whose unit cell contains two 
layers both of which can be made of isotropic chiral materials. The homogenized stack 
is a uniaxial bianisotropic medium (UBM) [16]. Numerical results on spectroscopic 
shifts are discussed in Section 3, followed by concluding remarks in Section 4. An 
exp(— iujt) time-dependence is implicit, with uj denoting the angular frequency. The 
free-space wavenumber, the free-space wavelength, and the speed of light in free space 
are denoted by k Q = Uy/e fi , \ = 2n/k , and c G = 1/ ^/fi e , respectively, with ji Q and 
s being the permeability and permittivity of free space. Vectors are in boldface, dyadics 
underlined twice; column vectors are in boldface and enclosed within square brackets, 
whereas matrixes are underlined twice and similarly bracketed. Cartesian unit vectors 
are identified as u x , u y and u z . The dyadics employed in the following sections can be 
treated as 3x3 matrixes [17]. 



2. Theory 

Let us consider a periodic multilayered stack whose unit cell is made of two layers of 
dissimilar materials labeled 1 and 2, as shown in Fig. 1. Both materials are isotropic 
chiral, with their Tellegen constitutive relations written in the frequency domain as [5]- 
[10] 

D = e Q e n E + ic'^nH 



B = /i G /i n H - ic Q ^ n E 



71=1,2, (1) 



where e n is the relative permittivity, fi n is the relative permeability, and is the chirality 
parameter. The thickness of the layer made of the n-th material is denoted by d n , so 
that 

fn = d n /{d 1 + d 2 ), 71 = 1,2, (2) 

is the volume fraction of the n-th material in the stack. 

Let the z axis of a Cartesian coordinate system be oriented normal to the layers. 
Provided that both layers in the unit cell are electrically thin, the periodic stack can 
be homogenized as a uniaxial bianisotropic medium (UBM) whose distinguished axis is 
the z axis [16]. The frequency-domain constitutive relations of the UBM are given as 

D = e [e±L+ 0|| - e±)u z u z ] «E 

B = Ho [n± I + (fx\\ - fi±)u z u z ] ' H 

where J is the identity dyadic, and there are six effective constitutive parameters: s\\, 
e±, fx\\, £||, and In the following subsections, we provide two different ways 
to determine these six scalars in terms of the constitutive parameters and the volume 
fractions of the two constituent materials in the unit cell. 



(3) 
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Figure 1. Schematic of the periodic multilayered stack of interest. 



2.1. 4 X 4-Matrix method 

Without loss of generality, wave propagation in the multilayered stack can be handled 
by using the spatial Fourier transform as 

E(x, y, z) = e(z, k) exp (ikx) dn 

H(x, y, z) = h(z, k) exp (ikx) dn 

whereby the fields are assumed not to vary along the y axis. The spatial frequency 
along the x axis is denoted by k. Substitution of this representation along with the 
constitutive relations (1) in the two Maxwell curl equations leads to the 4x4-matrix 
ordinary differential equation [18] 



(4) 



>M1 



[F{z,k) 



(5) 



in any layer made of the n-th material. Here, the column vector 

e x (z, k) 



[F(z, «)] 



e y (z, k) 
h x (z,K) 
h y (z, k) 



(6) 



represents components of the electric and magnetic fields that are tangential to the 
bimaterial interfaces in the multilayered stack, whereas the 4x4 matrix 



p(«)( 



K 











t^ot^n 


i^n/ Co 





l^ot^n 













i^n/ Co 
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+ 



K 



-i^nCo 



e n /fx 





"t^n/^-o 



— i^ n c Q 




n= 1,2, 



(7) 



contains all three constitutive parameters, the angular frequency, and the parameter k 
denoting propagation along the 

A similar exercise for wave propagation in the homogenized medium (i.e., the 
equivalent UBM) yields the 4x4- matrix ordinary differential equation [18] 

<l «)]=<[£**(«)] • [F(z,k)}, (8) 



dz 



where 



+ 



~i£±/c fi n± 

i£±/c Q -HoPx 

-e Q e± -i£±/c 

e Q e ± i£±/c 

-^||C -n\\/e 



e\\/fi -i£\\c 





(9) 



As both layers in the unit cell are electrically thin, we can invoke the long 
wavelength approximation [19] to set 



\E eqvt (K)} = fl 



+/ 2 



With 



A = {he 2 + / 2 ei)(/iAi 2 + / 2 /ii) - (A6 + Mi) 2 , 



(10) 



(11) 



this procedure yields the following expressions for the effective constitutive parameters 
of the stack: 



e± 
e \\ = 



A" 1 [/iei(e 2 Ai 2 - £D + f2e 2 (em - , 
: A" 1 [/i/ii(e 2 /i 2 - f|) + /2/^( £ i^i - fi)] , 

A- 1 [/i6(^ 2 - e 2 ) + / 2 6(^i - £ 2 )] • 



(12) 
(13) 
(14) 
(15) 
(16) 
(17) 



Thus, whereas the mixing of the constitutive parameters of the two layers of the 
unit cell is trivially simple in the transversely isotropic components (e±, fi±, of 
the constitutive dyadics of the equivalent UBM, the mixing is far richer in the axial 
components (e\\, fx\\, £y) of those dyadics. Our focus is on these axial components. 
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2.2. Method of boundary conditions 

A critical question that arises in any homogenization procedure is, whether the procedure 
preserves the boundary conditions across interfaces that are imposed by the Maxwell 
equations on the electromagnetic fields. In this subsection, we obtain the effective 
constitutive dyadics of the periodic multilayered stack by considering the boundary 
conditions on the electromagnetic fields. The obtained parameters are identical to the 
ones obtained by the procedure of Sec. 2.1, which is a matter of consistency. 

In the limit of very small layer thickness (compared to the wavelength of light), 
the electromagnetic fields can be assumed to be reasonably uniform across a unit cell of 
the periodic stack. But they have to satisfy the appropriate boundary conditions across 
the bimaterial interfaces. Suppose the fields in layer n, n e [1,2], of the unit cell are 
denoted by E n , H n , D„, and B n . Let us first consider the continuity of the tangential 
components (i.e., oriented parallel to the xy plane) of the E and H fields: 

Eltan = E 2 ton = (E tan ), (18) 
Ultari = H 2 tan = (H tan ), (19) 

where (Etan) and (H tan ) are the volume-averaged fields. Now considering the volume- 
averaged fields (D tan ) and (B tan ) defined as 

(D ton ) = flD Uan + f 2 ~D 2ta n = £o£±(Eton) + ic'^iHtan), (20) 

(Btan) = /lBi tan + f 2 B 2tan = flofl±(Htan) — ^C~^_L (E tan ) , (21) 

we obtain 

e± = fiei + f 2 e 2 , (22) 

A*± = ht*i + hto, (23) 

a = Mi + / 2 6 (24) 

for the equivalent UBM. These results are identical to those given by Eqs. (12), (14), 
and (16), respectively. 

Next, let us consider the continuity of the electromagnetic field components normal 
to the bimaterial interfaces (i.e., oriented along the z axis) as follows: 

(D norm ) ~Dlnorm D 2norm , (25) 

(B norm ) ~Q\ norm B 2norrrt , (26) 

(E norm ) fl^lnorrn ~\~ f2^-'2normi (27) 
(H-norm) /lHi norm + f 2 H. 2 norm- (28) 

Using the constitutive relations (3) of the homogenized material, we obtain the following 
three relations from the foregoing equations: 

_£i + f £ j_ 

ei^i - ii e 2 fi 2 - £ 2 

eifii - ii e 2 ^ 2 - & 



£|| 


£ m - 








e\\n\\ - 




en 









h^^ + h—^, (29) 
/iT-^ + Zs^h^ (30) 



h — ^ T2 +/ 2 — ^—5. (3i) 
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The solutions of these three equations are given by Eqs. (13), (15), and (17). Let us 
reiterate that, as all expressions derived in this section involve only the volume fractions 
and not the individual layer thicknesses, they are valid only in the limit of small layer 
thicknesses. 

The quantities e\\, and £y satisfy the relation 







£2 






Det 


in 


y>2 
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= 
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6 







This relation implies that both homogenization procedures yield only two of the three 
quantities e\\, and £y independently; the third can be obtained from Eq. (32). This 
feature of homogenization has been remarked upon earlier for dispersions of electrically 
small, isotropic chiral spheres in an isotropic achiral host material [20]. 



3. Numerical Results and Discussion 



In order to prove the premise of this paper, let us set material 1 to be achiral, i.e., £i = 0. 
Furthermore, without loss of generality, let both materials 1 and 2 be nonmagnetic in 
the Tellegen representation: ^ — /i 2 — 1-|| In that case, /i± — 1 whereas 

(/i£2 + jjgi) ~ Ml 

^R+MFJg- (33) 

Now, we can explore the effect of £ 2 on the sign of Re {e\\). 
Suppose further that material 1 is silver, so that [21] 

£l (A„) =5.7 + 0.4,- (^)V (34, 



where q = 1.6022 x 10~ 19 C is the charge of an electron and K = 1.0546 x 10~ 34 J s 
is the reduced Planck constant. This expression is valid for A D e [300, 900] nm. It has 
been used, for example, by Pendry [22] and also fits the data of Stahrenberg et al. [23] 
in the chosen spectral regime. For e 2 , we choose the Lorentz model 

£ 2 (A )=2 + —S!* — , (35) 

1 _ [ | _ »72A2 [ A2 I 

\\o) 2wc \X J 

wherein C 2 , X2, and 72 are constants. The values of these parameters chosen for the 
following studies are typical of solid materials; in addition, they keep the magnitudes of 
e 2 from becoming unphysically large in the resonance regime. 

Figure 2 contains plots of e\\ and £y with respect to free-space wavelength, when 
material 2 is nondispersive and nondissipative, and £ 2 = 0.05. The resonant frequency 
evident in these plots is due to the condition A = 0; hence, it is fixed by the volume 
fraction of silver (fi = 0.58 for this figure) and is affected by the dispersive properties 

|| An isotropic chiral material may be nonmagnetic in the Tellegen representation, but not in the 
Drudc-Born-Fedorov representation, and vice versa [20] . 
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of silver. The width of the resonance is proportional to the volume fraction of the 
dissipative material (i.e., silver). For physically realistic values of £2 10 -3 ), the 
resonance condition A = can be achieved only if the ratio Re(£2)/Re(£i) < 0. 

The resonance condition enhances £y. This becomes clear by noting that the 
maximum value of |£||| is many times larger than that of the chirality parameter £2 of 
the isotropic chiral constituent of the multilayered stack. In contrast, Eq. (16) indicates 
that £j_ is diluted in relation to £ 2 by the volume fraction of material 2. 

The resonance condition also enhances Im(£||) far above the imaginary part of 
the relative permittivity of silver. This is deleterious to wave propagation inside the 
multilayered stack. Still, we cannot help remarking that frequency regimes exist (e.g., 
A G ~ 500 nm) wherein £y is significantly enhanced whereas Im(£||) is significantly reduced 
in relation to Im(£i). 

Let us also note that the homogenized medium, a UBM, has the properties of 
a chiroplasma except for a gyrotropic term, because such a term is missing in both 
constituent materials [24]. It also presents an example of a homogeneous medium that 
displays chirality along with a negative real permittivity at optical frequencies, which 
has not been reported earlier, to our knowledge. 

Figure 3 presents a highly magnified view of e\\ in the resonance region (A ~ 
400 nm). The constitutive properties and geometric parameters are the same as in the 
previous figure, except that data is shown for both £2 = and £2 = 0.05. Clearly, 
the zero-crossing of Re(£ii) blueshifts by about 0.5 A, when £ 2 increases from to 0.05, 
furthermore, the same blueshift would occur even if £2 were to be replaced by —£2- This 
spectroscopic shift is measurable [15]. 

We have not shown spectral plots of /iy, because its magnitude is very close to unity. 
It does, however, evince a resonance, just like £y and £y and at the same frequency. 

Let us now move on to the situation where the relative permittivity of material 
2 is both dissipative and dispersive. In order to investigate the interplay of structural 
resonance (evident in Fig. 2) and intrinsic material resonance, we set A2 = 400 nm; 
furthermore, C2 = 10~ 3 and 72 = 10 13 rad s -1 . Figure 4 clearly shows the doubly- 
resonant characteristics of £11 and £m, arising from this interplay. Chirality enhancement 
is again evident in this figure, as also the enhancement of Im(£||). 

The double resonance in the spectra implies that there are three zero-crossings of 
Re(£||). Figure 5 presents a highly magnified view of e\\ near one of the three zero- 
crossings for both £2 = and £ 2 = 0.05. A blueshift of 0.25 A is evident, which, we 
reiterate, is measurable [15]. 

We also calculated the linear remittances (i.e., the reflectances and transmittances) 
of a 20-period multilayered stack that was about a wavelength thick. The presence of 
isotropic chirality resulted in the aforementioned blueshifts in the effective constitutive 
parameters of the stack, and manifested as small changes in the remittances at oblique 
incidence. Cross-polarized remittances on the order of 1% were seen when £ 2 7^ 0, but 
were of course absent for £ 2 = 0. At normal incidence, £2 played no role at all, which 
becomes clear from setting k = on the right side of Eq. (9); instead, the stack behaves 
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Figure 2. Real and imaginary parts of en, and real and imaginary parts of as 
functions of free-space wavelength, when fi = 0.58, C2 = 0, and £2 = 0.05. 
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Figure 3. Real and imaginary parts of en as functions of free-space wavelength, when 
fi = 0.58 and C2 = 0; £2 = (solid lines) and £2 = 0.05 (dashed lines). 



aS 9b dilute metal. 

Our focus on spectroscopic shifts at the boundary of the visible and the ultraviolet 
regimes is in accord with the fact that the chirality of many isotropic materials is evident 
in both regimes. Examples include suspensions of poly-L-glutamic acid [25], solutions 
of glucose [26] and several derivatives of cellulose [27], and chiral barbituric acid [28]. 
Flooding a multilayered stack containing void regions alternating with metallic layers 
with isotropic chiral fluids is a way to dynamically blueshift the zero-crossing of Re (e ( |). 

Let us note that we have used £2 = 0.05 in order to illustrate the spectral shifts due 
to chirality. Chirality paramaters of natural as well as synthetic materials available in the 
published literature are smaller. However, values of |^| ~ 0.1 at microwave frequencies 
have recently been obtained [14]. When such large chirality parameters would become 



Spectral shifts due to isotropic chirality 



10 



30 
20 

^=10 



300 350 400 450 500 550 

K (nm) 



600 



0.3 
0.2 


-0.1 
-0.2 



a; 



300 350 400 450 500 550 

X (nm) 



600 



50 
40 
3=30 
E 20 
10 




0.4 



0.3 



UU 1 



0.2 



0.1 



300 350 



400 450 500 

X (nm) 



550 600 



300 350 



400 450 500 

X (nm) 



550 600 



Figure 4. Real and imaginary parts of ey, and real and imaginary parts of 
as functions of free-space wavelength, when fi = 0.58, C2 = 10~ 3 , A 2 = 400 nm, 
72 = 10 13 rad s _1 , and £ 2 = 0.05. 
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Figure 5. Real and imaginary parts of eii as functions of free-space wavelength, when 
/i = 0.58, C 2 = 10" 3 , A 2 = 400 nm, and 72 = 10 13 rad s^ 1 ; 6 = (solid lines) and 
£2 = 0.05 (dashed lines). 



available at optical frequencies, even larger spectral shifts would be possible. Effectively, 
the shift of the resonance frequency can enable the use of such stacks as tunable high- 
pass filters for spectroscopy. The use of two such filters in tandem can enable a narrow- 
bandwidth source at almost any desired frequency. Dissipation is, however, a cause for 
deep concern. We are investigating the possibility of making infrared filters using layers 
of low-dissipation polaritonic crystals with negative real permittivity (such as LiTaC>3 
resonance at about 26.7 THz or silicon carbide in the mid-infrared regime) and isotropic 
chiral materials. Such filters could be effectively used for spectroscopy of molecular 
vibro-rotational levels. 
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To conclude, we derived the effective constitutive parameters of a period multilayered 
stack whose unit cell comprises a metallic layer and an isotropic chiral layer. We found 
that the stack behaves as an anisotropic metamaterial with resonant response that is 
determined primarily by the volume fraction of the metal. The presence of isotropic 
chirality can result in small shifts of the resonance frequency for reasonably large values 
of the chirality parameter; hence, the sign of one of the two effective permittivities of the 
stack can be switched. The spectral shifts in the dielectric properties can be potentially 
useful for spectroscopic purposes. 
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